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JEFFREY S. GERONIMO AND PLAMEN ILIEV 

Abstract. We give a complete characterization of the positive trigono- 
metric polynomials Q(0,ip) on the bi-circle, which can be factored as 
Q(6, <p) = |p(e lS , e lv )| 2 where p(z, w) is a polynomial nonzero for \z\ = 1 
and \w\ < 1. The conditions are in terms of recurrence coefficients 
associated with the polynomials in lexicographical and reverse lexico- 
graphical ordering orthogonal with respect to the weight A ^iq^ g ^ on 
the bi-circle. We use this result to describe how specific factorizations of 
weights on the bi-circle can be translated into identities relating the re- 
currence coefficients for the corresponding polynomials and vice versa. 
In particular, we characterize the Borel measures on the bi-circle for 
which the coefficients multiplying the reverse polynomials associated 
with the two operators: multiplication by z in lexicographical ordering 
and multiplication by w in reverse lexicographical ordering vanish af- 
ter a particular point. This can be considered as a spectral type result 
analogous to the characterization of the Bernstein-Szego measures on 
the unit circle. 



1. Introduction 

The factorization of positive polynomials as a sum of squares of polyno- 
mials or rational functions is an important problem in mathematics and led 
Hilbert to pose his 17th problem which was solved by Artin. In the case 
of trigonometric polynomials one of the simplest factorization results is the 
lemma of Fejer-Riesz which states that every positive trigonometric polyno- 
mial Q n {9) of degree n can be written as Q n {6) = \Pn(e ie )\ 2 where p n {z) is a 
polynomial of degree n in z. This result has been useful for the trigonometric 
moment problem, orthogonal polynomials, wavelets, and signal processing. 

Extensions of this result to the multivariable case cannot be generic as a 
simple degree of freedom calculation on the coefficients shows. Recently 
these results have been extended to two variable factorizations 

Qn,m(0,<P) = |Pn, m (e i0 ,e^)| 2 (1.1) 

where n and m are the degrees of Q n , m in 6 and ip, respectively, in the case 
when p ntm (z, w) is a polynomial of degree n in z and minw which is nonzero 
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for \z\ < 1 and \w\ < 1. This augments results obtained earlier by Kummert 
[I~8] (see also Ball [1]), Cole and Wermer [5], and Agler and McCarthy pQ 
(see also Knese [IE]). In particular, using the results of Knese [17] it is easy 
to see that, except for certain special cases, the polynomials p n ,m in <|l-lj) 
cannot be associated with the distinguished varieties defined by Agler and 
McCarthy [lj. Some extensions to more than two variables of the above 
results have also recently been obtain by Grinshpan et al [14], Bakonyi 
and Woerdeman [3], and Woerdeman [19] . In this paper we extend the 
results in [IT] in a different direction. We characterize completely positive 
trigonometric polynomials <3n,m(^;9 2 ) ; which can be factored as in (jl.ip 
where p n ,m(z,w) is a polynomial which is nonzero for \z\ = 1 and \w\ < 
1. The conditions can be written in a relatively simple form if we use 
the orthogonal polynomials in lexicographical and reverse lexicographical 
ordering introduced in [12] with respect to the weight 471 ^q 1 m ^ on the 
bi-circle. More precisely, in Theorem 12.41 we prove that (jl.ip holds if and 
only if certain matrices (which represent recurrence coefficients) K. n ^ m , K,\ m , 

V n ,m-, ^n,m satisfy the equations 

^m^iAJ^lmf = 0, for j = 0,1, . . . ,n- 1. (1.2) 

There are two important cases when equation (jl.2p holds: 

(i) The case when /C njm = characterizes the stable factorizations of 
Qn,m discussed in [TT] (i.e. equation (jl.ip holds with a polynomial 
Pn,m(z,w) which is nonzero for \z\ < 1 and < 1). 

(ii) The case when /C* m = characterizes the anti-stable factorizations 
of Qn,m- In this case (jl.ip holds with a polynomial p ntTn (z,w) such 
that z n p njm (l/z,w) ^ for \z\ < 1 and \w\ < 1). 

We derive several corollaries of the above result which are of independent 
interest. For instance, we characterize the Borel measures on the bi-circle 

for which the recurrence coefficients Ek,u ^k,l multiplying the reverse poly- 
nomials associated with the two operators: multiplication by z in lexico- 
graphical ordering and multiplication by w in reverse lexicographical vanish 
after a particular point, see Theorem 12.101 This can be considered as a 
spectral theory type result analogous to the characterization of Bernstein- 
Szego measures on the circle. We also show that in this case the space of 
orthogonal polynomials can be decomposed as an appropriate direct sum 
of two sets of orthogonal polynomials associated with the stable and the 
anti-stable factorizations described above, see Theorem 12.71 

The paper is organized as follows. In Section [2] we introduce the notation 
used throughout the paper including the recurrence formulas and state the 
main theorems. In Section [3] some preliminary results are proved and certain 
relations among the recurrence coefficients developed in [9] and their conse- 
quences are discussed. In Section [4] we prove the first main theorem which 
yields the factorizations (jl.ip withp„ )m (z, w) nonzero for \z\ = 1 and \w\ < 1. 
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In the forward direction, we use the Gohberg-Semencul formula, paramet- 
ric and matrix- valued orthogonal polynomials to show that if holds, 
then the recurrence coefficients i for the polynomials in lexicographical 
ordering associated with the weight ilT iq 1 ^ ^ vanish after a particular 

point. This leads to (|1.2j) . The heart of the proof in the opposite direction 
is based on a very subtle decomposition of the space of polynomials in the 
reverse lexicographical ordering as the sum of two subspaces possessing a lot 
of extra orthogonality properties. Using this decomposition, we construct an 
appropriate rotation on the space of polynomials which gives the polynomial 
Pn,m( z ,w) satisfying (jl.ip . All these constructions are missing in the stable 
case: the space decomposition is trivial (one of the subspaces is empty) and 
the rotation is simply the identity transformation. Thus, in our construc- 
tion, the polynomial p n ^ m {z } w) is no longer the first column of the inverse of 
the Toeplitz matrix associated with the trigonometric moments, but instead 
is a linear combination of the columns in the first block column of this ma- 
trix. One can use also the general theory of Helson and Lowdenslager |15| 
and the constructions in Delsarte et al [7] to obtain factorizations of positive 
functions Q(Q,cp) on the bi-circle. Note, however, that their approach works 
in a rather general setting and will provide (in general) non-polynomial fac- 
torizations of Q n ,m, even when (jl.ip holds with a polynomial p n ,m(z, w). In 
Section [5] we prove all remaining statements and corollaries. In Section [6] 
some examples are presented as illustrations of the main theorems. 

2. Statement of results 

2.1. Basic notations. We denote T = {z € C : \z\ = 1} the unit circle and 

T 2 = {(z,w) : \z\ = \w\ = 1}, 

the bi-circle (torus) in C 2 . Throughout the paper, we will use the parametriza- 
tion z = e l9 and w = e lLp , where 9, tp 6 [— 7T, it]. 

We consider moment matrices associated with the lexicographical order- 
ing which is defined by 

(k,£) <i ox (h,£i) o k < ki or (k = ki and £ < £\), 

and the reverse lexicographical ordering defined by 

(M) <rovlox (h,h) ^ (£,k) <i ex (£i,h). 

Both of these orderings are linear orders and in addition they satisfy 

(k,£) < (m,n) => (k + p, £ + q) < (m + p,n + q). 

Let n n,m denote the bivariate Laurent linear subspace sp&n{z k w l , —n< 
k < n, —m < I < m} and let £ be a linear functional defined on n n,m such 
that 

C{z- k w- 1 ) = cfc,, = C(z k w l ). 
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We will call c\. \ the (k, I) moment of C and C a moment functional. If 
we form the (n + l)(m + 1) x (n + l)(m + 1) matrix C njm for C in the 
lexicographical ordering then it has the special block Toeplitz form 



C n 



Co C_i 
C\ Co 

C n C n _i 



C_ n 
C-n+l 

Co 



(2.1) 



where each C/% is an (m + 1) x (m + 1) Toeplitz matrix as follows 



Cfc 



Cfcn c fc,-i 



c fc. 



Cfc,0 



-n, . . . , n. 



(2.2) 



Thus C„ im has a doubly Toeplitz structure. If the reverse lexicographical 
ordering is used in place of the lexicographical ordering we obtain another 
moment matrix C n ^ m where the roles of n and m are interchanged. We say 
that the moment functional C : IP'" 1 — > C is positive if 



C\p(z,w)p(l/z,l/w)} > 



(2.3) 



for every nonzero polynomial p(z, w) £ IT 
set p(z, w) 



flC[z,4 Here and later we 



p(z, w). It follows from a simple quadratic form argument 
that C is positive if and only if its moment matrix C n ^ m is positive definite. 
We now perform the Gram-Schmidt procedure on the monomials using the 
lexicographical ordering. The study of orthogonal polynomials on the bi- 
circle with this ordering was begun by Delsarte et al. [6] and extended in [12] . 
Given a positive definite linear functional C : n * — > C we perform the 
Gram-Schmidt procedure using the lexicographical ordering on the spaces 
span{z k w l : < k < n, < I < m} where n < N , m < N . Thus we define 
the orthonormal polynomials </>® m (z,w), 0<n<N,0<m< M, < s < 
m, by the equations 



-k„..-h 



£(4>n,m(z,w)z W 



0, < k < n and < I < m or k = n and < I < s, 



^(</ > n,m(^ ) ^)0n,m( 1 / 'z,l/w)) = 1, 



and 



E 

(fc,0<lcx(™,s) 



(2.4) 



(2.5) 



With the convention kn',m,s > 0, the above equations uniquely specify 4> s n m . 
Polynomials orthonormal with respect to C but using the reverse lexico- 
graphical ordering will be denoted by <p & nm . They are uniquely determined 
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by the above relations with the roles of n and m interchanged. Set 



Y6 m 1 




- z n w m - 


±m-\ 
K rn,m 


— K 


z n w m-l 


6° 




1 



(2.6) 



where the (m + 1) X (n + l)(m + 1) matrix K n ^ m is given by 



K 



k 



in,ra-l 
n,m—l 



n,m,m— 1 



iii,0 ,11-1,1(1 



,0,0 
. 0,0 



"'n.m.m— 1 



€ n,m,0 



(2.7) 



As indicated above denote 



~dj n ' 

Yn,m 




w m z n 








Yn,m 




w m z n-l 


— K 




dp 




1 


-Yn.m- 







(2. 



where the (n+1) x (n+l)(m+l) matrix K n ^ m is given similarly to (I2.7P with 
the roles of n and m interchanged. For the bivariate polynomials m (z, w) 



■f- 



above we define the reverse polynomials (j) s nm (z,w) by the relation 



Cm^w) = Z n W m 4j"{l/z,l/w). 



(2.9) 



With this definition d)^^ m {z,w) is again a polynomial in z and u>, and fur- 
thermore 



<f> njm (z,w) :-- 



m 

n,m 
m— 1 
n,m 







(2.10) 



An analogous procedure is used to define 



We use M m ' n to denote 



the space of all m x n matrices. In [12] it was shown: 
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Theorem 2.1. Giwen {$„ im } and {<& n ,m}, < n < N, < m < M, the 
following recurrence formulas hold: 

^n,m^n,m = Z*&n— l,m ~~ ^n,m^n—l,mJ (2.11a) 
^n,m + ^n,m-^«,«i(^n,m) ^ n,m = ^\i,m z ^n~l,mi (2.11b) 
I^m^r^m = ^n,m- 1 K-"n,m*&n—l,mi (2.11c) 

rl n,rn$n,m = W®n,m-l ~ ^n,m * Ll,m» (2. lid) 

^n.,m In,m^n,m ,m^n,m— 1j (2. lie) 

^ n,m = In,m^n,m + (r n )Tt ) ^n,m— 1' (2.11f) 

w/iere 

^n,m = (^n-l.m, J = ^n, m G M m+1 ' m+1 , (2.12a) 

A„, m = (z* n - 1)m , $ n>m ) G M m+1 ' m+1 , (2.12b) 

/C n ,m = (^.m-l.i-l,™) G M m '", (2.12c) 
T„, m = ($ n ,m-l, $n,m> G M m ' m+1 , (2.12d) 

/C^ = <«;<JWi, <^_ 1>m > G M™>", (2.12e) 
= (w^m-!, d> n , m ) G M m > m+1 , (2.12f) 

/n,m = (V» Km) G M m+1 >" +1 , (2.12g) 

= G M m+1 ' n+1 . (2.12h) 

Remark 2.2. From now on we adapt the following convention. For every 
statement (resp. formula) we will refer to the analogous statement (resp. 
formula) with the roles of z and w exchanged as the tilde analog. For 
instance, the tilde analog of formula (|2,lla|) is A niTn $ n>m = io$ njTO _i — 

^n,m * n ,m—X' 

Finally, we note that 

2.2. Main results. We say that a polynomial p(z,w) G C [z, to] is of de- 
gree (n, m) where n and m are the minimal nonnegative integers such that 
p(z, w) G n n,m . We say that the polynomial p(z,w) is stable if it does not 
vanish for \z\ < 1 and \w\ < 1. Similarly, for a trigonometric polynomial 
Q(9, ip) = p{e ld , e J</3 ), we define the degree as the ordered pair (n, m), where 
n and m are the minimal nonnegative integers such that p(z, w) G n n ' m . 
We can now state our main results. 

Theorem 2.3. For a positive moment functional C defined on the space 
jjn.m ^ e f n ow ,i n g conditions are equivalent: 
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(i) There exists a polynomial p(z,w) of degree at most (n, m), nonzero 
for \z\ = 1 and \w\ < 1, such that 

I r ik6 ilip 

£( z k w l )= / . — dOdtp. (2.13) 

v ' 4?r 2 7 |p(e* e ,e^)| 2 ^ v ; 

[ — 7T,7t] 2 

(ii) T/ie coefficients K n>m) K} nm , V n>m , f^ m satisfy 

^nA^lAJ(Pn,m) T = 0, for j = 0,1, ... ,n - 1. (2.14) 
Moreover, if the conditions above hold, we have 

\p(z,w)\ 2 =$ nym {z,w) T <& n , m (z,w) - $> n ,m-l (z, w) T $ n ^ m _i (z, w) 

=§ n m (z, w) T $ n ^ m (z, w) - $ n -l,m(z, w) T & n - 1<m (z, w), (2.15) 
for (z,w) £ T 2 . 

The polynomial p(z, w) in Theorem 12.31 can be computed from equation 
(|4,lip in Section HI which depends on the matrices U and V constructed 
from IC n>m , tC^ m , T n ^ m and r„ m in Lemma HTBI and Lemma HHJ see Remark 
14.51 for more details. 

As an immediate corollary of the above theorem and the maximum en- 
tropy principle [2] we obtain the first Fejer-Riesz factorization. 

Theorem 2.4 (Fejer-Riesz I). Suppose that Q(9,ip) is a strictly positive 
trigonometric polynomial of degree (n,m). Then Q{Q,ip) = \p(e l9 , e tLp )\ 2 
where p(z, w) is a polynomial of degree (n, m) such that p(z, w) ^ for 
\z\ = 1, \w\ < 1 if and only if the coefficients K n . m , K} nm , f n>m , t\ m 
associated with the measure ^JZq^q ^ on [— ft, 7 ?} 2 satisfy equation (|2.14p . 

Analogous results hold with the roles of z and w and n and m interchanged 
if the coefficients in the reverse lexicographical ordering satisfy the tilde 
analogs of equation (I2.14p (see equation (j2. 17b[) below). In the case when 
both sets of conditions hold we find: 

Theorem 2.5. For a positive moment functional C defined on the space 
n n ' m the following conditions are equivalent: 

(i) There exist stable polynomials p(z , w) andq{z,w) of degrees (n\, m\) 
and {n2,rn2) with n\ + n2 < n, m\ + r?i2 < m such that 

i r ike Hip 

C( z k w l ) = — / — — — d6dip. (2.16) 

[ — 7T,7r] 2 

(ii) The Coefficients K^n,m? ^n,m^ ^n,rri} ^n,m> ^-n,mj ^n,m Satisfy 
^nM^m^kmY ( )C l,m) T = °; for j = 0,1, ... ,U - 1, (2.17a) 

4, m [ri, m ri,J^i im = 0, for l = 0, 1,. . . ,m- 1. (2.17b) 
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As in Theorem 12.31 given the coefficients in the recurrence formulas, the 
polynomial p(z, w)z ri2 q(l/z, w) can be computed by (|4.1ip (see also Remark 
14. 5p . In view of equation (I2.16p . we say in the rest of the paper that a func- 
tional satisfying the equivalent conditions in the above theorem belongs to 
the splitting case. Theorem 12.51 can also be recast as a Fejer-Riesz factoriza- 
tion. 



Theorem 2.6 (Fejer-Riesz II). Suppose that Q(6,(p) is a strictly positive 
trigonometric polynomial of degree (n, m). Then Q(9, ip) = \p(e ld ,e lip )q{e~ t9 ,t 
where p(z,w) and q{z,w) are stable polynomials of degrees (n-i,mi) and 
(tt-2, 772-2) respectively, with n\ + n-2 = n, m\ + mi = m if and only if the co- 
efficients /Cn,m; ^n,m; ^n,m, ^n, m ; ^n,m, ^nm associated with the measure 
iJ 9 Q{e, v ) on t -71 "' 71 "] 2 satis fv equations (|2.17p. 

In this case when the equivalent conditions in Theorem 12.51 hold we have 
the following structural theorem. 

Theorem 2.7. Suppose that (|2.16[) holds, where p(z,w) and q{z,w) are sta- 
ble polynomials of degrees (ni,m\) and (ri2, ra^), respectively. Let <& p kl {z,w) 
and Q Q kl (z,w) be the the (vector) polynomials orthogonal with respect to the 



measures 



dd dp 

47r 2 |p(e i; «,e^)| 2 



and 



dd dip 



*p{z,w) 



i-K 2 \q(e w ,e^)\ 2 



, respectively. Then 



(z,w) 



K 2 ,m 2 -l( Z > W ) 



(2.18) 

Moreover, if we set n = n\ + ni and m = mi + rri2, then there exist unitary 
matrices U € M m ' m , V £ M m+1 > m+1 such that 

- z n2 q(l/z,w)^ p numi ^(z,w) 

V( Z ) w)w T 



n,m—l 



(z,w) 



(2.19a) 



and 



p(z, w)z n2 q(l/ z, w) 
z<*q(l/z,w)3? ^w) 



lt(z,w)w m2 ^ m2 _ 1 (z,l/w)_ 

Roughly speaking, the above theorem allows to decompose the space of 
orthogonal polynomials associated with the functional in (I2.16P as a sum of 
the two extreme cases: 

• the stable case when q(z,w) = 1; 

• the anti-stable case when p(z,w) = 1. 

As a corollary of the the proof we obtain also the following characterizations 
of these situations. 

Corollary 2.8. For a positive moment functional C defined on the space 
pjn,m ^ e j \\ ow i n g statements hold. 



(2.19b) 
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(i) There exists a stable polynomial p(z, w) of degree at most (n, m) such 
that 

C( z k w l ) = — / ————dOd<p (2.20) 
v ; 4vr 2 J \p(e te , e^)\ 2 Y y ' 

[— w,tt] 2 

if and only iffC njjn = 0. Moreover, we can take p(z,w) = (f>™ m (z,w) . 

(ii) There exists a stable polynomial q(z, w) of degree at most (n, m) such 
that 

i f ik9 ilip 

C( z k w l ) = -^ / , , .„ . ^ dOdtp (2.21) 
v ; 4vr 2 J \q(e- te ,e^)\ 2 r 1 ; 

[-7r,7r] 2 

if and only if K} n m = 0. 

As a consequence of the above corollary, we obtain a simple character- 
ization of the functional which are tensor products of functionals on the 
circle. 

Corollary 2.9. Let C be a positive moment functional on the space n n ' m . 
Then, there exist a positive functional C z defined on span{z fe : \k\ < n} and 
a positive functional C w defined on span{w/ : |/| < m} such that C(z k w l ) = 
C z (z k )C w (w l ) if and only if K n<m = K\ m = 0. In this case, <j>™ m {z,w) = 
a{z)j3{w), where a(z) and /3(w) are stable polynomials of degrees at most n 
and m, respectively and 

1 r pike ilip 

C( z k w l ) = —, / , . * . d6d<p. (2.22) 
v ; 4vr 2 J \a(e ie )P{e^)\ 2 * 1 ; 

[-7T,7r] 2 

Finally, the above results can be used to completely characterize the mea- 
sures on T 2 for which the corresponding coefficients E^i and E k: i vanish after 
a particular point. 

Theorem 2.10. Let [i be a positive Borel measure supported on the bi-circle. 
Then fi is absolutely continuous with respect to Lebesgue measure with 

, dO dip 

M 4vr 2 |p(e ie \e^)q(e- i( > ',e^)| 2 ' 1 ' 

where p(z,w) and q(z,w) are stable polynomials of degrees (ni,mi) and 
(ri2, vn-i), respectively, with n\ + < n, mi + m2 < m if and only if 

E kjl = and E k j = for all k > n + 1, l>m + l. (2.24) 
Moreover, in this case we have 

Ek,i = 0, for k > n + 1, I > m — 1 and E^j = 0, for k > n — 1, I > m + 1. 

(2.25) 
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3. Preliminary results 

3.1. Connection between bivariate and matrix orthogonal polyno- 
mials. The vector polynomial <& n>m (z, w) denned in (12.6P can be written as 



®n,m(z,w) 



w 



w 

m—l 



1 



(3.1a) 



where (z) is a unique (m + 1) x (m + 1) matrix polynomial of degree 
n in z. Similarly, the vector polynomial & ntm (z, w) defined in (|2.8p can be 
written as 



z' u 



1 



(3.1b) 



where ^^(w) is a unique (n + 1) x (n + 1) matrix polynomial of degree m 
in w. The recurrence relation (j2.11a|) and its tilde- analog are equivalent to 
the recurrence relations for the matrix- valued polynomials {&™(z)} n >o and 

{$mH}m>0- 

We will also need the following Christoffel-Darboux formula, which is a 
tilde analog of formula (4.1a)-(4.1c) in [12] 



tr tr , tr tr 

$n,m(2, ^) $n,m (^1,^1 ) T ~ $ n-l,m(z, w) $ n— l,m 

(Zl, Wl) T 
~l,m(^)' w ) ^n— l,m(^l > W\) 



(3.2) 



= (1 - rowi)$ nim (z, w) T $ nim (zi, wi), 
and its corollary (see equation (4.2) in [12]) 



,m(^-l)'^l) ^n,™,— 1 ,m—l (zi,wx) 



(3.3) 



3.2. Relations among the coemcients. We list below different relations 
among the coefficients defined in (|2.12p needed in the paper. 

The tilde analog of formula (3.52) on page 811 in [12] can be written as 
follows 

ri+i./ffc+i.j =f fc>l f i,i + ih,iEk+i,i(il,i) T ^ ^ 
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We also need formulas (3.1), (3.4) and (3.6) from [9]: 

T^Ek+^Tl,) 1 +/C fc ,,(/Ci I ) J J (3.5a) 

Fk,iEk+i,iIk,i = Ak+i,i-i]Ck+i,l ~ fck,lFk,i> (3.5b) 

iljEwpijf = (rt+i,i) TA l"+i,i-i - ?lM,i) T - ( 3 - 5c ) 

Recall that if E^i = then Aki = is the identity + x (Z + l) matrix. 
Using this fact and the above formulas, we see that the following lemma 
holds. 

Lemma 3.1. If 

then 



Ek+l,l = and E k+1 ^ = 0, 


(3.6) 


&k,l(fik,l) T = °) 


(3.7a) 


pi pt _ pt pi 
1 k+1,1 1 k+l,l ~ 1 k,l L kp 


(3.7b) 


JCk+1,1 = fck,lTk,h 


(3.7c) 


(fcl+l,l) T = rl,l(fck,lf ■ 


(3.7d) 



3.3. Stability criterion. Throughout the paper we will use several times 
the following fact: a polynomial p(z,w) is stable (i.e. non-vanishing for 
\z\ < 1 and \w\ < 1) if and only if 

• p(z,w) ^ for \z\ = 1 and \w\ < 1, and 

• p(z,w) 7^ for \z\ < 1 and \w\ = 1. 

The above criterion is a simple corollary from the well-known stability cri- 
teria for bivariate polynomials, see for instance |S]. 

4. One sided stable polynomials 
In this section we prove Theorem! 



4.1. Proof of the implication (i)=^(ii) in Theorem 12.31 Assume first 
that the conditions in Theorem 12. 3( i) hold, i.e. the moment functional C is 
defined on n n,m by 

1 f e ik8 e iltp 



£(zV) = 4^ J w f ^*W dB * p ' 

[— 7T,7r] 2 

where p(z,w) is of degree (re, m) nonzero for \z\ = 1 and < 1. (4.1) 



We can use (|4.ip to extend the functional C on the space of all Laurent poly- 
nomials C[z, z -1 , w, w~ 1 ]. Thus we can define vector polynomials i(z, w) 
for all k,l G No- 



12 



J. GERONIMO AND P. ILIEV 



For every fixed z = e %9 E T, we denote by L the corresponding positive 
moment functional on the space C[tu,u7 _1 ] given by 



\p(e ie ,e^)\ 2 



dip. 



(4.2) 



Similarly, for a polynomial (j)(z, w) of degree (k, I) we can fix z = e %9 on the 
unit circle and consider the corresponding polynomial 4>(e l9 , w) of degree I 

in w which depends on the parameter 9. We will denote by <fi (e , w) the 
reverse polynomial of <p(e l9 ,w), i.e. we set 



c/) e (e i8 ,w) =w l 



'AM- 



Lemma 4.1. Suppose that (|4.ip holds. Then with respect to C 9 we have 



p(e l0 ,w) ±{w l :l> 0}, 
p^{e l9 ,w) ±{w l :l< to}, 



and \\p(e l9 , w)\\ ■ 
Proof. We have 



p e {e tB , w) 



1. 



(4.3a) 
(4.3b) 



C 9 {w l p{e i9 ,w^j 



1 

2k 



e ilip p(e i9 ,e il f) 
\p(e i9 ,e^)\ 2 



dip 



l-i 



2tt J p(e l9 ,w) 



dw = 0, 



for I > by Cauchy's residue theorem, establishing (|4.3ap . The second 
orthogonality follows bv_a similar computation. The assertion about the 

norms of p(e l9 ,w) and p 9 (e t9 ,w) is straightforward. □ 



We would like to construct now polynomials {<pf (w)}i>q orthonormal with 
respect to C . From Lemma 14.11 it follows that we can take 



= w l ~ m p°(e w ,w), for I > to. 



4- 



(4.4) 



Let us denote by C 9 the (I + 1) X (I + 1) Toeplitz matrix associated with 
C 9 , i.e. if we put c 9 = L (w~i) then 







r 9 1 

c -l 




C l C 


■ c 9 




c 9 c 9 

H c i-i 


r 9 
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Recall that we can use the coefficients of the orthonormal polynomial <f>f(w) 
to compute the inverse of C?-, via the Gohberg-Semencul formula [131 The- 
orem 6.2, page 88]. Explicitly, if we set 



3=0 



(4.5) 



then 



1-1/ 



6 



i-i 



'l-l 



o 
o' 



'l-l 



J J 



2-1 



o 



(4.6) 



Lemma 4.2. Suppose that (|4.1|) holds for all (k,l) G Z 2 . Then 

E k l = for k>n + \ and I > m - 1. (4.7) 

Proof. Note that for fixed I > m — 1, the matrix polynomials {^jU^OlfcXJ 
defined in Subsection 13.11 are orthonormal on [— tt, n] with respect to the 
matrix weight ^=Cf, i.e. 

± J\ l k (e ie )Cf[^(e ie )]U9 = 5 kj I l+1 . 

From the theory of matrix-valued orthogonal polynomials it will follow that 
E k l = for k > n + 1 if we can show that (Cf)^ 1 is a (matrix) trigonometric 
polynomial in 9 of degree at most n. This follows immediately from (j4.4j) . 
(14311 and M. □ 



(4.8) 



Lemma 4.3. Suppose that equation (j4.7|) holds. Then 



K 



k.i 



f i ft 



(JC\ ; ) T = 0, /or all j > 0, fe > n, Z > m. 



Proof. From Lemma 13.11 we see that equations (13. 7j) hold as long as k > n 
and I >m. First, we would like to show by induction on j £ No that 



^k+j,l(^k+j,l) 



k.l 



f 1 ft 



(/C^) T , for k>n, l> m. (4.9) 



If j = 0, the above statement is obvious. Suppose now that (14. 9ft holds for 
some j > 0. From (|37TB|) it follows that 

U'+l 



f l ft 



f 1 



r 



k+l,l 



J ~ 



t 
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Using the above formula we find 

i+i 



K 



f i ft 
1 fc,r fc/ 



kl) 



A:./ 



'-a ft 
/c+i,r fc+i,z 



/C 



fc+ij 



r 



i ft 

k+l,l L k+l,l 



{fc\ + n) T (by equations ()3.7c[) and (|3.7dp ) 

= fck+i+j,i(K-k+i+j,i) T (by the induction hypothesis), 

establishing (|4.9p for j + 1 and completing the induction. From (|3,7ap we 
see that the left-hand side of (|4.9p is equal to leading to (|4.8p . □ 



Proof of the implication (i)=t>(U) in Theorem \2. 31 The proof follows imme- 
diately from Lemma 14.21 and Lemma 14.31 □ 



4.2. Proof of the implication (ii)=^(i) in Theorem 12.31 The key in- 
gredient of the proof in the opposite direction, which also explains the con- 
struction of the polynomial p(z,w), is the following lemma. 

Lemma 4.4. Let £ be a positive moment functional defined on n n ' m . Sup- 
pose that there exist unitary matrices U 6 M n,n and V 6 M n+ 1 > n+1 such 
that 

>(1) 



and 



V*$n t m(z, 



*«-i, m (*,uO 



S*n-l,m(*» w ) 
*i-l,m(^™) 



(4.10a) 



(4.10b) 



where ^^li (z,w) is an nj -dimensional vector whose components are poly- 
nomials of degrees at most (n — 1, m) with n\ + nz = n, and m (z, w) is a 
polynomial of degree at most (n,m). Then 

<t -= 



p(z, w) = il> l m (z, w) = z n w m ^ m (l/z, 1/w) (4.11) 

is a polynomial of degree at most (n,m), nonzero for \z\ = 1, \w\ < 1 and 
equations ()2. 13j) and (|2,15p hold. 

Proof. From equations (|4.10p and (14. lip it follows that 

$„ im (z,-u;)$ nim (2;i,ii;i) T - $ n _i im (z, to) $ ra— l,m 



for zzi = 1, 



(4.12a) 



and 



(z, w) T $ njm (zi,u;i) - 

= p (z,U>)^p (zi,t/Ji) 



$n-l,m(^,w) $n-l,m(zi, Wl) 
for 2ffi = 1, 



(4.12b) 
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where p (z, w) 
obtain 



z n w m p(l/z,l/w). Plugging equations f)4. 12|) in f|3.2|) we 



p(z, w)p{l/ z,w\) — wwip(z, w) p (1/z, Wl) 

= (1 - WWl)$ n>rn (z, w) T $ nm (l/z, Wl). 



(4.13) 



Using the last equation we can prove that p(z,w) is nonzero for \z\ = 1 
and | it; | < 1. Recall first that the vector polynomials ^ n ,m(z, w) can be 
connected to the matrix polynomials <&™(z) via (13. lap . Moreover the matrix- 
valued orthogonal polynomials {^k,m( z )}k=o constructed in Subsection 13.11 
are orthonormal with respect to the matrix inner product 



(A,B) = C(AM m (w)B^), 
where M m {w) is the (m + 1) x (m + 1) Toeplitz matrix 



(4.14a) 



M m (w) 



w 



w 

m—1 



W 



w 



W 
1 



W 



w'" 

m—1 



W 



-m+1 



1 



(4.14b) 



In particular, from the theory of matrix-valued orthogonal polynomials we 
know that det[3>™(z)] / for \z\ > 1. This implies that 



(z,w) is a nonzero vector for \z\ = 1 and w£C. 



(4.15) 



Suppose first that p{zq,Wq) = for some \zq\ = 1 and |ioo| < 1. Then using 
()4. 13|) with z = zq and w = w\ = Wq we obtain 

-\w \ 2 \p(z , W )\ 2 = (1 - \w \ 2 )$ n>m (zo, W ) T ® n , m (z , W ). 

Since the left-hand side of the above equation is < and the right-hand side 
is > 0, we see that $ njm (zo ; ^o) must be the zero vector, which contradicts 

(gup . 

Suppose now that p(zo,wo) = for some \zq\ = 1 and \wq\ = 1. Then 
p(zo,Wo) = and therefore equation (14.131) with z = zq, w = wo and 
w\ / wq gives 



®n,m(zO,Wo) ^n.m^O^l) = for all Wl ± W , 

which implies that $ n ,m(zo,wo) is the zero vector leading to a contradiction, 
thus proving the required stability for p(z,w). 

Note that equation (I2.15[) follows easily from ()4.12a|) and ()3.3|) . Thus, 
it remains to prove that equation (|2.13j) holds. Let us denote by pi(z) the 
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coefficient of w l in p(z,w), i.e. we set 

m 

p(z,w) = ^2pi(z)w l . 

1=0 

Then a straightforward computation shows that for \z\ = 1 we have 
p(z, w)p(z, w\ ) — wwi*p(z, w) i p(z, Wl) 



(4.16) 



1 — WW\ 



[1 



w ■ ■ ■ w 







/ r Po (z) 



\ 



o 



.Pm(z) ••• p (z). 



Po(z) Pl(z) ... p rn {z) 



O 



Po(z) 



o 



Pm(Z) 



Pl(z) ••• p m {z) 

Prom (|3.1ap we see that 



Pm(z) 
O 



Pi{z) 



Pm{Z) 




1 

Wl 

wT 



(4.17) 



<S> n , m (z,w) T $ n , m (z,wi) = [1 w ••• w m ]J m ^(z) T ^(z)J, 



1 

Wi 



UK 



(4.18) 



where J m = [Si tTn -j]o<ij<m- From equations (|4.13p . (|4.17p and (|4.18p it 
follows that for \z\ = 1 we have 



Po{z) 
Pi(z) 



O 



Po{z) Pi(z) ... p m {z) 



.Pm(z) ■■■ Po(z)] I O 

O' 



Pm{Z) 



pi(z) ■■■ 
J m ^{z) T ^{z),l 



Pm{z) 



Po{z) 

p m {z) ••• Pl{z) 



O 



Pm{Z) 





(4.19) 



Since p(z, w) is nonzero for \z\ = 1 and \w\ < 1 we see that for fixed z = e 
on the unit circle, (^(w) = w m p(e td ,1/w) is an orthonormal polynomial 
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of degree m with respect to the (parametric) moment functional C 9 , with 
moments 



1 



4 = /:>-') = i- 



-ilip 



2ir J \p(e i9 ,e i f)\ 2 



dip, for |Z| < m. 



(4.20) 



From Gohberg-Semencul formula (see [131 Theorem 6.1, page 86]) it follows 
that the left-hand side of equation (14.19P is the inverse of the Toeplitz matrix 



rib 
<-'.» 



C C_i 



-m+1 



Since = I m+1 and J m C 6 m J n 



-m-l ' ' ' c 

(Cfj) T , equation (|4.19p gives 



[$™(z) f <S>n(zT\ where z = j e . (4.21) 

From the theory of matrix- valued orthogonal polynomials we know that the 
matrix weight on the right-hand side of (14.211) generates the same matrix- 
valued orthonormal polynomials Q<k<n an d therefore 



7T 

C(z k M m (w)) = i- / e lke C e n 



d6 for — n < k < n. 



From the first row of the last matrix equation we find that for — n < k < n 
we have 



1 

2^ 



jke. 



,d8. 



(4.22) 



The proof of ([2Tl3|) follows at once from equations (j4T20]) and (|4T22|) . □ 



Remark 4.5. To complete the proof of Theorem 12.31 we need to show that 
equation ()2.14j) implies the existence of unitary matrices U and V such that 
equations ()4.10p hold. For moment functionals satisfying ()2. 16j) the existence 
of such matrices follows easily from the tilde analog of formulas (j2. 19|) . In 
general (for one-sided stability) the construction of U and V is the content 
of the next two lemmas. 

Note also that if we know V we can compute explicitly p(z, w) in equation 
(|2.13p from (|4.10bp and (|4.1ip . In Theorem 12. 31 we gave the simplest formula 
for \p(z, w)\ 2 , which involves only the orthogonal polynomials. However, one 
can easily extract from the proof of Lemma [4.41 other formulas which can be 
used in practice to compute p(z, w). For instance, we can use (|4.12a|) (which 
is stronger than (|2.15p ). or setting w\ = in (|4.13p we obtain 

p(z,w)p{l/z,0) = $„, m (z,JK) T i^(l/z,0), 

which gives p(z, w) up to a factor depending only on z. 
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Lemma 4.6. Let JC and K, 1 be m x n matrices, and let r = rank(/C) ; r 1 = 
rank(/C 1 ). Then the following conditions are equivalent: 

(i) /C(/C 1 f = 0; 

(ii) We have 

JC = USU\ K> = U l S l U T , (4.23) 

where U, U 1 £ M m,m , U 6 M n,n are unitary and S, S 1 are m x n 
"diagonal" matrices with block structures of the form 



(4.24a) 





si 




s = 






















and 



S 1 



s 1 

S r i 



(4.24b) 



with positive s\, . . . s r , s\, . . . , s \ and r + r 1 < n. 
Remark. Note the the condition r + r 1 < n implies S'(S' 1 ) T = 0. 

Proof. We focus on the implication (i) =>- (ii), since the other direction is 
obvious. Consider yl = lOK, and -B = (K}) T K, 1 . Note that A and B are 
hermitian n x n matrices such that AB = BA = 0. Hence, there exists an 
orthonormal basis (iix,U2, • • • , u n ) for C n which diagonalizes A and B, i.e. 



K)K,Uj = XjUj, {K}) T K, l Uj = /XjUj and 



(4.25) 



Let U be the unitary matrix with columns u\, . . . ,u n . From ([4.25P we see 
that 



Xj = \\fCuj\\ > and 

Moreover, the fact IC{K}) T 
that 



lO-uAr > 0. 



Ai, A2, • • • A r are positive and A r+ i = • • • = A 
/i n _ r i + l, . . . ,/J-n are positive and \i\ = /X2 = • 



implies that XjfXj = for all j. Suppose now 

0; 



Mr, 



0. 



Set 



Si ■ 

1 



'A, 



J j V r"n—r L +j 

Consider the sets of vectors 
K,Ui 



for i 
for j 



1,2, 
1,2, 



■it,- 



: % 



1,2, 



(4.26a) 
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and 



T-i 



j+n—m 



J 



m 



r l + l, 



, m 



(4.26b) 



^ j—m+r 1 ) 

Using (I4.25P it is easy to see that T and T\ are orthonormal sets of vectors. 
Extending the set T to an orthonormal basis for C m and constructing a 
matrix with columns these vectors we obtain a unitary m x m matrix U ; 
extending the set T\ to an orthonormal basis for C m and constructing a 
matrix with columns these vectors we obtain a unitary m x m matrix U . 
With these matrices one can check that (|4.23p holds. □ 

Remark 4.7. If know that ICiK 1 ) 7 = and K)K} = (which are satisfied 
by the matrices JC = /C„ im and /C 1 = /C^ m in Theorem 12. 5h then we can 
choose U = U l in equation (|4.23p . Indeed, using the notations in the proof 
of Lemma 14.61 we see that vectors in T are perpendicular to the vectors in 
T\ . Extending the orthonormal set T UT\ to an orthonormal basis for C m 
we can construct a unitary matrix U = U 1 with columns these vectors. 



Lemma 4.8. Let G be an n x n matrix. Suppose that 

(G k ) id = for alii = 1,2,... , r, j = n - r 1 + 1, n - r l + 2, . 
and k = 1, 2, . . . ,n — 1, 



,n, 



(4.27) 

where r + r 1 < n. Then, there exists a unitary n x n block matrix £ of the 
form 





r I r 








£ = 





* 













I r l 



such that the matrix £^G£ has the following block structure 

&G£ 



* 


" 


* 


* 



(4.28) 



(4.29) 



where the zero block in equation ()4.29p above is an n± x n<i matrix with 
n\ > r, n2 > r 1 and n\ + ri2 = n. 

Proof. Let {ei, e2, . . . , e n } be the standard basis for C n , and let 
Wo = span{e n _ r i +1 ,e n _ r i +2 , . . . ,e n }. 

Consider the space 

= W + GW + • • • + G n_1 W - (4.30) 

By Cayley-Hamilton theorem, W is the minimal subspace of C n which is 
G-invariant and contains Wq. From equation ()4.27p it follows that 



spanjei, e2, . . . , e r } C W . 



(4.31) 



Let 
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• {vi,V2, ■ ■ ■ ,v ni } be an orthonormal basis for W 1 - which extends 
{ei,e 2 , . . . ,e r }, i.e. vj = ej for j = 1,2, ... ,r; 

• {v ni+ i, v ni+ 2, • • • , v n } be an orthonormal basis for W which extends 
{e n _ r i +1 ,e n _ r i +2 , ■ ■ . ,e n }, i.e. = ej for j > n- r 1 . 

Then the unitary matrix £ with columns v\ , V2 , • • • , v n will have the block 
structure given in (j4.28[) , and the G-invariance of W implies equation (14.29|) . 

□ 

Proof of the implication ( ii J=^> (i) in Theorem \2.S\ Applying Lemma l4.6l with 
/C = /Cfj, yyi and /C = m 

we see that there exist unitary matrices U, U 1 £ M m,m , U £ M n,n such that 
equations (|4.23p . (|4.24j) hold. Moreover, applying Lemma 14.81 with 

G = U^t 1 ft U 

we see that U can be modified (if necessary) , so that equations (14.231) , (|4.24p 
hold and 



G = fj+f 1 ft U 

n.m n.m^ 







(4.32) 



where the zero block in the equation above is an n\ x ri2 matrix with n\ > 
r = rank(/C nim ), ri2 > r 1 = rank(/C* m ) and n\ + n 2 = n. 

Replacing /C n>m and m in the tilde analogs of equations (|2.11cp and 
(|2.11d|) with the expressions given in (|4.23j) we find 

t^tf n : m^n,m 

= C/t$ n _ l m _ 5 T j7t$ n)TO _ lj (|4.33b ) 

and 



r/tf 1 $ 

u 1 n,m^n,m 



With ?ii and n 2 fixed above, we will use the following notation: for an 
ra-dimensional vector ^, we denote by (resp. ^j>( 2 )) the vector which 
consists of the first m (resp. the last n 2 ) entries of the vector Thus if we 
set 

*n-l,m = U^n-^m, (4.34) 

then the vector m can be represented in the block form 



n— l.m 



(1) ' 
n— l.m 
(2) 

n— l.m 



(4.35) 



With this choice of a unitary matrix U, we want to show that there exists a 
unitary matrix V such that (|4.10b|) holds. Since the bottom n 2 rows of the 
matrix S T are equal to 0, we see from equations (|4.33b ) and (|4.35j) that 



(t/tf r 



^n, m ) (2) = ^ n \ m . 635t) 

Similarly, since the first ri\ rows of the matrix {S l ) T are equal to 0, we see 
from equations (|4.33d ) and (|4.35p that 



(2) 



(^1* m * n , m )« = Z* 



(1) 

n— l.m' 



d436H) 
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Equations fj4.36j) show that the entries of the vector polynomials z^n-i m 

~ (2) 

and ^ n _i m are linear combinations of the entries of the vector polynomial 
& n ^m and therefore, they are orthogonal with respect to C to all monomials of 
degree at most (n, m— 1). Moreover, since U is unitary, it follows from (|4.34[) 

and (14.35P that the entries of each of the vectors z^^-i m ano ^ *n-l m f° rm 
orthonormal sets of polynomials of degrees at most (n, m) with respect to C. 
Finally, from equations (|4.32p and (|4.36p we see that the entries of the vector 

z ^n-i m are perpendicular to the entries of the vector ^^-l m- Therefore, 
all the entries in the vectors z^S^-i m an< ^ *n-i m f° rm an orthonormal set of 
n polynomials, which can be extended by adding a polynomial rfi™ m (z, w) to 
an orthonormal set of polynomials of degree at most (n, m), perpendicular to 
all polynomials of degree at most (ra, m — 1). The transition matrix between 
this set and the orthonormal polynomials {<p s n m (z, w)} s =o,i,...,n is a unitary 
matrix whose transpose is a unitary matrix V satisfying equation (I4.10b|) . 
The proof now follows from Lemma l4.4i □ 

5. Proofs of the theorems in the splitting case 

5.1. Proof of Theorem 12.51 The proof of implication (i)=^(ii) follows eas- 
ily from Theorem 12.31 and its tilde analog. Indeed, note that if fl2. 161) holds, 
then we have 

1 f e ik8 e iltp 

£( z k w l )= / ——^ — ^--rrdOdip (5.1a) 

[— 7T,7r] 2 

1 f pikd ilip 

dOdcp (5.1b) 



(5.2a) 



(5.2b) 



4vr 2 J \Q(e id ,e^)\ 2 

[— 7T,7r] 2 

where 

P(z,w) = p(z,w)z n2 q(l/z,w) is a polynomial of degree at 
most (n, m), stable for \z\ = 1 and |u;| < 1, 

and 

Q(z,w) = p(z,w)w m2 q(z, 1/w) is a polynomial of degree at 
most (n, m), stable for \z\ < 1 and \w\ = 1. 

Therefore, equation (|2.17ap follows from Theorem l2.3l (i)=^(ii) for the poly- 
nomial P(z, w) and equation (|2.17bp follows from the tilde analog of The- 
orem [23] (i)=^(ii) for the polynomial Q(z,w). Conversely, suppose that 
equations (|2.17p hold. Then, by Theorem 12.31 (ii)=^(i) and its tilde analog, 
we deduce that there exist polynomials P(z,w) and Q(z,w) of degrees at 
most (n, m) such that equations f)5. 1 1) hold and 

(a) P(z,w) is stable for \z\ = 1 and \w\ < 1 

(b) Q{z,w) is stable for \z\ < 1 and \w\ = 1. 
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Without any restrictions, we can assume that P(z,w) and Q(z,w) are not 
divisible by z and to. From equation (12.151) for P(z,w) and Q(z,w) we see 
that \P(z, w)\ 2 = \Q(z, w)\ 2 for all (z, w) £ T 2 which implies that in the ring 
of Laurent polynomials C[z, z^ 1 , w, w^ 1 } we have 

P{z, w)P{l/z, l/w) = Q{z, w)Q(l/z, l/w). (5.3) 

Suppose now that we factor P(z, w) into a product of irreducible factors 
Pj(z,w) in C[z,w]. Likewise we can factor Q(z,w) into a product of irre- 
ducible factors qi(z,w) in C[z, w\. Using (15. 3p we see that we can factor 
P(z,w)P(l/z,l/w) = Q(z, w)Q(l/z, l/w) in C^z" 1 ,^,^" 1 ] in two ways 
as a product of irreducible factors 

P(z,w)P(l/z, l/w) = Y[ Pj (z, w)p j (l/z, l/w) 

TT 1 ~ (5 ' 4) 

= [[qi(z,w)qi(l/z,l/w) = Q{z,w)Q(l/z,l/w). 
l 

Since C[z, z , w, w^ 1 ] is a unique factorization domain, it follows that for 
every j, there exists a unique I such that exactly one of the following holds: 

(I) pj(z,w) and qi(z,w) are associates in <C[z, z , w, it; -1 ]; 
(II) pj(z,w) and q~i(l/z,l/w) are associates in C[z, z -1 , w, to -1 ]. 

Note that the units in C[z, z~ l , w, to -1 ] are of the form cz s w r , where c ^ 0, 
s,r £ Z. Thus, we see that if (I) holds then with the normalization chosen 
above (P and Q are not divisible by z and w) we must have pj(z,w) = 
cqi(z,w) where c is a nonzero constant. From properties (a) and (b) of the 
polynomials P(z,w) and Q(z,w) we deduce that pj(z,w) = cqi(z,w) ^ 
when \z\ = 1, |to| < 1, and likewise pj(z,w) = cqi(z,w) / when \z\ < 1, 
| to | = 1. This shows that if (I) holds, then pj(z,w) ^ when \z\ < 1 and 
\w\ < 1. 

If (II) holds then pj(z,w) = cz s ->w r i q~i(l/z, l/w) where c / 0, and Sj,rj 
are the minimal nonnegative integers for which z s ^w Tj qi(l/z,l/w) belongs 
to C[z, w]. It is easy to see that pj(z,w) and qi(z,w) have the same de- 
gree (sj,rj). From property (a) of the polynomial P(z,w) we deduce that 
z Sj Pj(l/z,w) = cw r J q~i(z,l/w) t^O when \z\ = 1, \w\ < 1. From property (b) 
of the polynomial Q(z, w) we conclude that z s ipj(l/z, w) = cw r ^qi(z, l/w) ^ 
when \z\ < 1, |u>| = 1. Thus we see that if (II) holds, then the polynomial 
z Sj Pj(l/z,w) has no zeros when \z\ < 1 and |tt;| < 1. 

Let J\ (resp. J2) denote the set of indices j for which (I) (resp. (II)) holds. 
Then the polynomials p(z, w) = IX/eJi Pj( z i w ) andq(,z, w) = Y\j € j 2 z s ipj(l/z,w) 
satisfy the conditions in Theorem I2.5f i). completing the proof. □ 

5.2. Proof of Theorem 12. 71 For the proof of Theorem 12.71 we summarize 
first some basic properties of the vector orthogonal polynomials &k,i(z,w) 
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associated with a moment functional C of the form 

■ d9 dip, 



\p(e ie ,e^)\ 2 



where p(z,w) is of degree (re, m) nonzero for \z\ < 1 and \w\ < 1. (5.5) 

We define as usual ^(z^w) = z n w m p(l/z, l/w). 
Lemma 5.1. If (|5.5p holds then 

C(p(z,w)z' k w' 1 ) = for all k E Z, Z > 0, (5.6a) 

= ° /or a// fc<n, IgZ. (5.6b) 



Proof. The proof of (|5.6a|) follows immediately by computing first the w 
integral and by using equation (|4.3ap in Lemma 14.11 The proof of (|5.6b|) 
follows by a similar computation, by evaluating first the z integral. □ 

Lemma 5.2. Suppose that (|5,5p holds. Then the vector polynomial ^ n ,m(z, w) 
has the following block structure 



V(^,w) 

$n,m-l(z,w) 



(5.7) 



Moreover, 

C{$ n ,m-i{z, w)z~ k w' [ ) = for allk<n, l> 0. (5.8) 
Proof. Equation (15.71) follows from Theorem 7.2 in [12]. Plugging (15. 7p and 



its tilde analog in (13. 2p we obtain the following identity 

p(z,w)p(zi,w 1 ) - < p(z,w) < p(z 1 ,w 1 ) 



= (1 - WWi)^ n . m -i(z,w) ® n ,m-l(zi,Wl) 
+ (1 - ZZ\) $ 

n— l,m 

(z,w) ^n-i^^^u;^ 

Thus, if take z = z\ on the unit circle the last term above will vanish and 
we can rewrite the equation as follows 

p(z, w)p(l/z, w~i) — ^p~(z, w)z~ n w™p(z, 1/tDi) 



1 — WWl 



l(z, w) T $ n!m -i(z,Wi). 



Using the matrix- valued polynomial <J>™ 1 (z) defined in (|3.1a.j) and its re- 
verse 
by 



m— 1 / 



z n <&™ (l/z) T we can replace in the last equation <& n)m _i{z, w\) 



z- n %™- l ( 



w 



-m— 2 
1 



and therefore we obtain 

p(z, w)z n p(l/ z, idi) — V( z > w)w~i l p(z, l/wi] 

1 — WWl 



$ nim _ 1 (^, U ;) T ¥^- 1 (z) 



\T 



<~ 2 



(5.9) 
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Let us denote by S(z, w, w±) the function on the left-hand side above. Note 
that we can rewrite S(z, w, wi) as follows 



_ A(z, 1/w, wi) <_ B(z,l/w,wi) 
b{z, w,w\) = p{z, w) h p (z, w) — 



it- 



ic 



(5.10) 



where 



and 



A(z, l/w,wi) 



B(z, 1/w,w-l) 



z n p{l/z, wi) — z n p(l/z, 1/w) 

1/w — W\ 

w~ m p(z, w) — w™p{z, 1/wi) 



1/w — Wl 

are polynomials in z, 1/w and w\ of degrees at most n, m — 1 and m — 1, 
respectively. Thus, there exist 1 x m vectors A m (z,l/w) and B m (z,l/w) 
whose entries are polynomials in z and 1/w of degrees at most n and m—1, 
respectively, such that 



A(z, 1/w, wi) = A m (z, 1/w) 



-m-l 

h 

-m-2 



and B(z, 1/w, wi) = B m (z, 1/w) 



-m-V 
-m-2 



Combining the last equation with equations (|5.9p and (|5.10p we see that 



A m (z,l/w) B m (z,l/w) 
p{z,w) p (z,w) = <E> 



w 



in 



^^(z^ft^Hzf. 



From the theory of matrix-valued orthogonal polynomials we know that 

det(<I>™ _1 (z)) 7^0for|2:| < 1. Therefore, the entries of the matrix [<!>™ _1 (-z) T ]~ 1 
are analytic functions on closed unit disk \z\ < 1 and we have 



r I A m (z,l/w) <_ B m (z,l/w) 

$n,m-l{z,W) =\P{z,w) p{z,w) 

1 w w 



m-l/ \T 



z) J 



Equation (|5.8|) follows immediately from the last equation and Lemma 15.11 

□ 



Proof of Theorem 2.7. The block structure of the vector polynomials given 



in equation (12.18P follows immediately from equation (15. 7ft in Lemma 15.21 
Let us denote by C p and C q the positive moment functionals corresponding 
to the stable polynomials p(z,w) and q{z,w), i.e. 



and 



C p (z k w l ) 



£ q (z k w l ) 



1 f e ike e ily> 
4^2 J \p(e iS ,e il fi)\ 2 

[ — 7T,7r] 2 



d9 dtp 



1 

4^ 



\q(e i9 ,e i! fi)\ 2 



d9 dip. 



To prove that there exists a unitary matrix U such that ()2.19ap holds, it is 
enough to show two things: 
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(i) the entries of the vector polynomial on the right-hand side of (|2.19ap 
form an orthonormal set of polynomials of degrees at most (n, m— 1) 
with respect to £; 

(ii) the entries of the vector polynomial on the right-hand side of (I2.19ap 
are orthogonal with respect to C to all polynomials of degrees at most 
(n — 1, m — 1). 

Clearly, the entries of the vector polynomial on the right-hand side of (|2.19ap 
are polynomials of degrees at most (n, m — 1) and it is easy to see that they 
all have norm 1. The fact that they are mutually orthogonal follows from 
the following direct computation 

£(z- n2 q(z, lM^^l/z, l/w) <p(z, w)w^- l ^ m2 _ 1 {z, l/w) T ) 



1 

4^2 



dz 

z 



p(z, w)z U2 q(l/ z, w) 

T LT 

= o, 

since the w-integral is zero by Cauchy's residue theorem. Thus, it remains 
to check (ii). Below we compute the inner products with the monomials 
z k w l where 0<k<n— l = ni + n2 — 1 and < I < m — 1 = m\ + m2 — 1. 
For the first mi entries on the right-hand side of (|2.19ap we obtain 

C(z^q(l/z,w)^ nitmi ^(z,w)z~ k w- 1 ) 

1 f Kx^-M™)^ 2 ^™' 1 dz dw 



4tt 2 J \p(z,w)\ 2 q(z,l/w) z w 

T 2 

/ z n 2 -k \ 

= £ P [K^-iM l/w) ) = o, 

by equation (|5.8p in Lemma [5.21 and the computation for the last m-i entries 
on the right-hand side of ()2.19ap is similar. The fact that there exists a 
unitary matrix V such that (|2.19bp holds can be established along the same 
lines. □ 

Remark 5.3. The decomposition in Theorem 12.71 can be naturally con- 
nected to a decomposition of a Christoffel-Darboux type formula. Indeed, for 
a polynomial h(z, w) let us consider the corresponding Christoffel-Darboux 
kernel 



T ht ^ h(z,w)h(l/z,rj) - h(z,w)h(l/z,r)) 
L {z,w;rj) = = . 

1 — W7] 

Using the notations in Theorem 12 .7\ we set q w (z,w) = z n2 q(l/z, w) and 
h(z,w) = p(z,w)q w (z,w). Then it is easy to see that 

L h (z, w; rj) = q w (z, w)q w (l/z, 7])L p (z, w; r]) + ^(z, rf)L q (z, w; if). 

The point now is that the polynomials p(z, w) and q(z, w) are stable for 
\ z \ < 1) | if | < 1 and therefore the corresponding kernels possess a great 
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many orthogonality relations (see for instance [10]) which can be used to 
prove equations (|2,19p and thus give an alternate proof of the implication 
(i)=>(ii) in Theorem 12.51 

It is a challenging problem to find a direct algebro-geometric proof of 
the implication (i)=>(ii) in Theorem 12.31 (i.e. if we have stability only with 
respect to one of the variables). If we use the notations in Theorem 12.31 
and Lemma 14.41 the heart of the problem is the following: start with a 
polynomial p(z,w) which is stable for \z\ = 1 and \w\ < 1 and give an ex- 
plicit description (or prove existence) of the spaces Hi and H2, where Hj 

is the space spanned by the entries of the vector polynomials m ( z i w ) 
in Lemma 14.41 These spaces must be mutually orthogonal and must satisfy 
additional extra orthogonality properties in view of equation (I4,10bj) . Since 
the orthogonality relations can be expressed in terms of residues, construct- 
ing bases for these spaces amounts to an interesting interpolation problem 
on a zero-dimensional variety, which involves appropriate zeros of p(z, w) 
and *p(z,w). Equivalently, this would give a subtle decomposition of the 
Christoffel-Darboux kernel associated with p(z,w). 

5.3. Proofs of Corollaries [2j5] and f279l 



Proof of Corollary 1 2. 81 The statement in (i) is proved in [12\ Theorem 7.2] 
but we sketch it briefly below since it follows easily from the constructions 
in this paper. If (|2.20p holds then the defining relation (|2,12cp for /C n)m 
and Lemma 15.21 show that tC n ^ m = 0. Conversely, suppose that JC n ^ m = 0. 
Equation (|2.11cp and its tilde analog imply that 



4>n,m( z , w ) 
$n,m-l(z,w) 



, $ njm (z,w) 



J n,m 



$n-l,m(z,w)_ 



(5.11) 



Using the second equation above and Lemma 14.41 (with U and V being 
the identity matrices), we see that equation (|2.20p holds where p(z,w) = 
(j) m m (z,w) is stable for \z\ = 1 and \w\ < 1. Since cj)™ m (z,w) = (fr™ m (z,w) 
we can use the first equation in (|5.1ip and the tilde analog of Lemma 14.41 to 
deduce that p(z,w) is stable also for \z\ < 1 and \w\ = 1, which shows that 
p(z,w) is stable for \z\ < 1 and \w\ < 1 completing the proof of (i). 



Suppose now that (|2.2ip holds. Applying Theorem 12. 71 and its tilde analog 
we see that there exist unitary matrices U £ M m ' m , U £ M n,n such that 

$ n , m -l(z, W) = W^UW^T^Z, l/w), $ n _l, m (z, w) = Z^U ^ n _^ m {l / 'z , 

Plugging these formulas in the definition ()2.12ep of K\ m we find 

Kim = U (K,m-lWz, ^ ^n-l,m(lA, ™)) ■ 
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Note that the inner product in the expression above gives the matrix }C n , 



for the measure 



d9 dip 



, and therefore is zero from the first part of the 

0. 



47r 2 |cj(e 

corollary. Thus, K\ m 

Conversely, suppose now that K} n m = 0. From (|2.11dj) and its tilde analog 
we see that there exist unitary matrices V G M m+l,m+l , V G M n+l,n+1 such 
that 



From Lemma [4.4l we deduce that (|2.2ip holds with q(z, w) 



z n j. n 
41 Y n,; 



(5.12) 

(l/z,w), 



which is stable for \z\ = 1 and \w\ < 1. We want to show next that ip™ m {z, w) 
and ij}™ m (z,w) are equal up to a unimodular constant, i.e. 



where 



1. 



(5.13) 



Note that if we can prove the above equation, we can use the tilde analog 
of Lemma 14.41 to deduce that q(z,w) = ez n tp^ m (l/ z, w) is stable for \z\ < 1 
and | it; | = 1, thus proving that q(z,w) is stable for \z\ < 1 and < 1. 

The proof of (|5.13p follows from the characteristic properties of ip™ m {z, w) 
and ip™m( z i w )- Indeed, from the first equation in (|5.12p it is easy to see 
that ^™ m {z,w) is the unique (up to a unimodular constant) orthonormal 
vector in n™' m such that 

V>™ m (z, w) _L {z k w l : < k < n - 1, < / < m} U {z n w l :l<l<m}. 

Similarly, from the second equation in (|5.12p we see that t/)™ m (z,w) is the 
unique (up to a unimodular constant) orthonormal vector in n n,m such that 

i>™ m (z, w) J_ {z k w l :0<k<n, < Z < m - 1} U {z k w m : 1 < k < n}. 



n,m 



z,w) and ip™ m (z, w) establish 

□ 



The above characteristic properties of tp. 
(|5.13p . thus completing the proof. 

Proof of Corollarv \2.9l Assume first that C(z k w l ) = C z (z k )C w (w l ). If we 
denote by {ctk( z )}o<k<n the (one- variable) polynomials orthonormal with 
respect to C z and by {A( w )}o<z<m the (one- variable) polynomials orthonor- 
mal with respect to C w , then it easy to see that 



<5>n,m{z,w) = a n (z) 



p m -x{w) 



L AM 



and $ nim (z,w) = j3 m (w) 



' a n (z) 
a n -i(z) 

. Oi (z) . 



From these explicit formulas and the defining relations ()2. 12c|) . (|2.12ep for 
K. n ^ m and K} n m it easy to see that /C njTn = /C„ m = 0. 

Conversely, suppose that /C n>m = K\ m = 0. Note that if h(z, w) is a 
polynomial of degree (k, I) such that h(z, w) and z k h(l/z, w) are stable, then 
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h(z,w) is independent of z (i.e. k = 0). Using this observation, Corollary 
12. 81 and arguments similar to the ones we used in the proof of the implication 
(ii)=^(i) in Theorem 12.51 we see that <f> ™ m (z, w) = a(z)/3(w), where a(z) 
and j3(w) are stable polynomials of degrees at most n and m, respectively 
and that equation (12.221) holds, completing the proof. □ 

5.4. Proof of Theorem 12.101 With the measure dp, we will associate the 
positive moment functional C defined on C[z, z~ 1 , w, t/j -1 ] by 

C{z k w l ) = j z k w l dn. (5.14) 

T 2 

First suppose that equation f|2.23|) holds where p(z, w) and q(z, w) are 
stable polynomials of degrees (ni,mi) and (772,777,2), respectively, with m + 
ri2 < n, mi + 777,2 < m. Then we can represent C as in equations (|5,ip where 
P(z,w) and Q(z,w) are given in equations (|5.2p . Using Lemma HT21 and its 
tilde analog, we see that equation (|2.25p holds. To complete the proof of 
the theorem, it remains to show that equation (|2.24p implies the existence 
of stable polynomials p(z,w) and q(z,w) of degrees (ni,mi) and (ri2,m2), 
with n\ + n2 < n, mi + 771,2 < m such that (I2.23P holds. From Lemma 14.31 
and its tilde analog we see that 

(/C^) T = 0, for all j > 0, k>n + 2, l>m + 2, 

(5.15a) 

K\ tl = 0, for all j > 0, k > n + 2, Z > m + 2. 

(5.15b) 

By Theorem 12.51 we deduce that there exist stable polynomials p(z,w) and 
q(z, w) of degrees (m, mi) and (77-2, 7772), with n\ + 712 < 77 + 2, mi + 7772 < 
777 + 2 such that (12TT61) holds for all (k, I) satisfying \k\ < n + 2, |/| < m + 2. 
Moreover, from equations (I2.15P we see that 

p(z, w)p(l/z, l/w)q(l/z, w)q(z, 1/w) 

= &n+2, m +2(z,w) T $n+2,m+2(l/z,l/w) - $„ + i >m+2 (z, w) T ^ n+lt m + 2 (l/z, l/w) 
= $ n +2,m+2(z, w) T '$ n+2 , m +2(l/ 'z, l/w) - § n +2,m+l{z, w) T $ n+ 2,m+l (l/z, l/w). 

(5.16) 

Recall that if E^ i = then A^j = and therefore by (I2.11ap we obtain 

$ k ,l{z,w) = z$ k - lt i(z,w). (5.17a) 

Similarly, if E/. ; = then 



K 



k.l 



f 1 ft 
1 k,l L fej 



(5.17b) 
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Using (I2.24|) we see that equations (|5.17p hold for all k > n+l and I > ra+1, 
which combined with (I5.16P shows that 

p(z, w)p(l/z, l/w)q(l/z, w)q(z, l/w) 

= $ k>l (z, w) T $ k j(l/z, l/w) - wf&k-uil/z, l/w) (5.18a) 

for k > n + 2, I > m, and 

p(z, w)p(l/z, l/w)q(l/z, w)q(z, l/w) 

= $ M (z, wf^l/z, l/w) - $ k ,i-x(z, wf^T^l/z, l/w) (5.18b) 

for k > n, I > m + 2. From Theorem 12.31 we see that equation (12.161) holds 
for all k, I € Z which establishes (|2.23[) . It remains to show now that in fact 
n\ + ?i2 < n and mi + m<i < m. To see this, we will use the following two 
observations: 

(i) If p(z,w) and q(z,w) are stable polynomials of degrees (ni,mi) 
and (712,777,2), then P(z,w) = p(z,w)z n2 q(l/z,w) is a polynomial 
of degree (774 + 722,777-1 + 7772) which is not divisible by z and w (i.e. 
P(0,tjj) ^ and P(z,0) 0). 

(ii) If P(z,w) is a polynomial of degree (no,777-o), which is not divisible 
by z and w such that P(z,w)P(l/z, l/w) £ n n ' m , then hq < n and 
777o < m - 

From (i) we see that the polynomial P(z,w) = p(z,w)z n2 q(l/z,w) is a 
polynomial of degree (ni + 772, 7771 + 7772) which is not divisible by z and ttj. 
From equation (I5.18ap with k = n + 2, / = 777 and equation (I5.18bj) with 
k = n, I = m + 2 we see that 

p(z, w)P(i/z, i/w) g n n+2 ' m n n n ' m+2 = n n ' m , 

which combined with (ii) completes the proof. □ 



6. Examples 

We now consider some examples that exhibit the properties of the theo- 
rems proved earlier. 

6.1. One-sided. Our first example will be a polynomial of degree (2,2) that 
is stable for \z\ = 1, \w\ < 1. We will construct the polynomial using the 
algorithm given in [12]. Setting n ,o = 1, ti 2 ,o = V 4 , ""-1,2 = j^jii u 2,2 = 

- ^iT a2) , u-2,2 = with (-3 + v / 13)/2 < a < (3 + >/l3)/2 and 

= 0, for (i,j) € {(0,1), (1,0), (0,2), (-1,1), (1,1), (1,2), (-2,1), (2,1)} 
we construct the orthogonal polynomials up to level (2,2). In this case we 
find using Maple or Mathematica that 

- 



2V15(1 
15(1 + a 2 ) 





-1 
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,2 



^2,2 



2,2 



15(1 - a 2 ) 
60a 

1 




1 2 




V3c 
3(l+a 2 ) 



r 1 

1 2,2 



2,2 



y^l-a 4 ) 
120^ 





2^3d 
3Vc 

1 



n Vc 8y / 15(l-a 2 ) 2 
15(l+a 2 )v7 







15(l+a 2 ) 



and 



2,2 



'Bd 




15 5 



(l-a 4 )(l-a 2 )v^ = 

4a/ 30a VJ 

5{l+a 2 )^c 4y / l5(l-a 2 ) 2 
2a^ 150^// 



where c = 14a 2 — a 4 — 1, d = 11a" — a 



4 - 1, d = 11c 2 <•' 

2 i 1 1„4 a„& i 7„2 



1, e = 55 + 190a 2 + 55a 4 , 

/ = 11 + 38a 2 + 11a 4 , g = 4a 4 + 7a 2 + 4. 

It is not difficult to see that £2,2(^2 2) T = ® = ^2 2^2 2> an< ^ usm S Maple 
or Mathematica we find 



while 



^2,2? 2,2? 2,2(^2,2) T = 
(1-a 2 ) 2 



iff pi ft jrl 

H2 i 2,2i 2 ,2^2,2- 15fl(1 + a 2 



^0. 



(6.1) 



Using Remark 14.51 we find a candidate for p(z,w) is 

p(z,u>) =(4a(l - a 2 )w 2 - 3(1 + a 2 ) 2 )z 2 + 3((1 - a 4 )u> 2 + 3a(l + a 2 ))z 
- 13a(l - a 2 )w 2 + 12a 2 . 

Using the Schur-Cohn test it is not difficult to see that p(z, w) is nonzero for 
\z\ = 1 and I it; I < 1. Applying Lemma [4.61 (see also Remark 14.71) we see that 
equations (Q3D and (04]) hold with r = r 1 = 1, Sl = J^Sj , »i = j^, 



C/ 1 = C/ 



1 

1 



and U = —= 



2 -1 
-1 -2 



Next we find ^1,2(2, w) in (|4.10ap by computing 



10a 



Aazw 2 —w 2 — a 2 w 2 + z — a 2 z 
-2azw 2 - 2w 2 - 2a 2 w 2 + 2z- 2a 2 2 



We look for a unitary matrix V such that equation (I4.10b[) holds with n\ 
n 2 = 1. 
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This uniquely specifies V (except the first column, which can be multiplied 
by an arbitrary complex number of modulus 1) as 



V 



\/3a 









V~c 


2V3(l+a 2 )Vd 


6a(l+a 2 ) 








v7 


Vd 




2 v / 3(l+a 2 ) 


v7 


v7 


v7 



The first entry of the vector polynomial ^2,2(25 w) = V^$2,2( 



Z.W) IS 



30a Vd 



P (z, 



Thus, Lemma H3 shows that p(z,w) = ~^~qP( z ' w )- 



6.2. Splitting case. The second example we will consider illustrates The- 
orem |23J In this case we chose uo,o = 1, u~i,i = a, ui,i = b, u 2i o = ab = uo <2 
and Ui d = 0, e {(0, 1), (1, 0); (1, 2), (-1, 2), (2, 1)', (-2, 1), '(2, 2), (-2, 2)] 
where —1 < a < 1 and — 1 < b < 1. Using the algorithm given in |12| we 
find that 

'a 



2,2 



r 1 

1 2,2 



^2,2 = 
^2,2 = 

o vT 
o 





0' 

b 







1 



r 



l o o' 

o VT^¥ o 



2,2, 



f 1 
1 2,2- 



It is easy to check that equations (|2.17p are satisfied. Moreover, we find that 
P(z,w)P(l/z,l/w) = ^ 2 {z,w) T \ 2 {l/z,l/w)-^ 2 {z,w) T \ 2 {l/z,l/w) 

where 

P(z, w 



2) 



(1 — bzw)(z — aw) 



^(l_ a 2 )(1 _ 6 2) 

is stable for \z\ = 1 and |u>| < 1. It is easy to see that the polynomial P(z, w) 
above is the unique polynomial (up to a multiplicative constant of modulus 
1) of degree at most (2, 2) which is stable for \z\ = 1 and \w\ < 1 and which 
satisfies (|6.2p . Finally, note that 

P(z,w) = p(z,w)zq(l/ z,w) 

where 



p(z,w) 



1 — bzw 



and q(z,w) 



1 



azw 



are stable polynomials. We can obtain all this also by following the steps of 
Example 1. Indeed, we see that we can take U = U = I 2 the identity 2x2 
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matrix and 

"0 1 0' 
V= 1 
1 

The first entry of the vector polynomial ^2,2 ( z , w) = ^'$2,2(2, w) is P(z, w). 
Note that if a = then ^2,2 = and the functional is in the stable case, 
while if b = then /C2 2 = and the functional is in the anti-stable case. 
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